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Diffusive transport of light in three-dimensional disordered Voronoi structures 
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The origin of diffusive transport of light in dry foams is still under debate. In this paper, we 
consider the random walks of photons as they are reflected or transmitted by liquid films according 
to the rules of ray optics. The foams are approximately modeled by three-dimensional Voronoi 
tessellations with varying degree of disorder. We study two cases: a constant intensity reflectance 
and the reflectance of thin fllms. Especially in the second case, we find that in the experimentally 
important regime for the film thicknesses, the transport-mean-free path l* does not significantly 
depend on the topological and geometrical disorder of the Voronoi foams including the periodic 
Kelvin foam. This may indicate that the detailed structure of foams is not crucial for understanding 
the diffusive transport of light. Furthermore, our theoretical values for l* fall in the same range as 
the experimental values observed in dry foams. One can therefore argue that liquid films contribute 
substantially to the diffusive transport of light in dry foams. 

PACS numbers: 82.70.Rr, 42.25.Dd, 42. 68. Ay, 05.60.-k 
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I. INTRODUCTION 

The interaction of light with matter is a highly inter- 
esting subject with many facets. For example, multiply 
scattered light in turbid media is a very complex topic. 
However, it becomes treatable when, after a sufficiently 
large amount of scattering events, light reaches its dif- 
fusive limit [l|. Diffusing light is used for biomedical 
imaging f2| and in difFusing-wave spectroscopy (DWS) it 
is able to monitor dynamic processes in turbid materials 
[31 . Recently, even clear evidence for strong localization 
of light was given [J|. In colloidal suspensions and ne- 
matic liquid crystals [5| the respective scattering events 
and diffusion of light are well understood. But why are 
aqueous foams [flHl: as used, e.g, in shaving cream, tur- 
bid? 

Aqueous foams consist of gas bubbles separated by liq- 
uid films p, Q- In dry foams, where most of the liquid 
has drained out, the bubbles are deformed to polyhedra. 
The so-called Plateau rules state that always three films 
meet in the Plateau borders, which form a network of liq- 
uid channels throughout the foam, and that four of these 
borders meet at tetrahedral vertices. Foams are visibly 
opaque although both the gas and liquid components are 
transparent. 

Precise light-scattering experiments show that light 
transport reaches its diffusive limit in foams |8l4lOl| , which 
means that photons perform a random walk. However, 
the mechanisms underlying this random walk are not well 
understood. One suggestion is that scattering from the 
Plateau borders is responsible for the diffusing light _9]. 
The transport-mean- free path /*, over which the pho- 
ton direction becomes randomized, is then predicted as 



I* ex Hj \fe, where H is the average bubble diameter and 
£ is the liquid volume fraction. However, data rule this 
out in favor of the empirical law [^, [llj 
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The authors of Ref. [llj state that this may imply sig- 
nificant contributions from scattering from vertices [l2| 
or films [13l-[l7|. Novel transport effects, such as total 
internal reflection of photons inside the Plateau borders, 
are studied in [TllfTsj. 

In recent years we have explored how liquid films in 
combination with ray optics determine light transport 
in dry foams [l3l - [l7j . In a step- by-step approach, we 
employed three model foams: the honeycomb structure, 
two-dimensional Voronoi tessellations 19ll. and finally the 
three-dimensional Kelvin structure [7|, |20[ . Apparently, 
the 'periodic Kelvin structure is highly idealistic. In his 
momentous experimental study of bubble shapes, Matzke 
did not find even a single Kelvin cell [21|. Therefore, in 
this paper, we extend our studies towards real dry foams. 
We introduce topological and geometrical disorder based 
on a three-dimensional Voronoi foam model to investigate 
the infiuence of disorder. Note that a Voronoi foam does 
not obey Plateau's rules, and indeed all films of a Voronoi 
foam are flat. Thus an accurate sample foam provided by 
the Surface Evolver software [22| clearly better describes 
the geometry of a real foam. Nevertheless, we deliber- 
ately concentrate on the Voronoi tessellations with fiat 
films in order to greatly simplify our simulations of light 
transport, which are still very time consuming. More- 
over, a comparison of the results presented in this paper 
with future theoretical studies of diffusive light transport 
in exact models of real foams, will highlight the role of 
both their geometry and their curved films. 

Cells in a foam are much larger than the wavelength of 
light, thus one can employ ray optics and follow a light 



TABLE I: Topological and geometrical characteristics of the cells in the Kelvin foam, our Voronoi foams 1-5, real foams [25 
simulated foams 26], and the Poisson Voronoi tessellation (PVT) [13|. The symbols are explained in the text. 



Kelvin Foam Foam 1 Foam 2 Foam 3 Foam 4 Foam 5 Real Foam Simulation PVT 

0.050197 0.126543 0.164922 0.312005 0.486383 

0.00684 0.01768 0.02340 0.04330 0.07433 

0.011179 0.028264 0.037067 0.069619 0.122512 

0.752663 0.748894 0.745840 0.727823 0.686616 

14.0 14.0 14.0 14.007812 14.188232 

5.142857 5.142857 5.142857 5.143335 5.154229 

0.0 0.0 0.0 0.028259 0.930633 

0.979593 0.979593 0.979593 0.979593 0.981509 1.028764 0.365-0.42 
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beam or photon as it is reflected by the liquid films with a 
probability r called the intensity reflectance. We perform 
extensive simulations to determine the diffusive limit of 
light for two models. In model I, we choose a constant 
intensity reflectance r to explore the effect of disorder. 
The essential result of our Monte- Carlo simulations is 
summarized in the empirical formula for the transport- 
mean- free path Z*, 



II. MODEL FOAM 



Voronoi 



{r)^0.63H- — -{l-bi+b2r). (2) 



The main behavior of 'voronoi(^) ^^ governed by the factor 
(1 — r)/r. Quite remarkably, this factor is also found in 
the honeycomb, two-dimensional disordered Voronoi, and 
three-dimensional regular Kelvin structures. However, 
there is a small but systematic deviation from (1 — r)/r, 
described by the last factor in Eq. ([2]) with < 6i < 0.1 
and 62 ~ 0.1. Both constants show a slight dependence 
on disorder in the Voronoi foam. In model II, we use 
the intensity reflectance of thin films with its significant 
dependence on the incident angle and on d/A, where d is 
film thickness and A the wavelength of light. Our theo- 
retical values for the transport-mean-free path I* fall in 
the same range as the experimental values for the driest 
foams in Ref. [9]. This shows that liquid films are im- 
portant for the understanding of photon diffusion in dry 
foams otherwise we should have obtained much larger 
transport-mean-free paths. We also observe that topo- 
logical and geometrical disorder do not change the qual- 
itative behavior oil* . Even quantitative changes are not 
very pronounced. Thus, our extensive numerical simu- 
lations support a conclusion that we already drew from 
our previous two dimensional investigations |14| : the de- 
tailed structure of a Voronoi foam is not important for 
understanding the diffusive limit of light transport. 

Our article is organized as follows. In Section|lT]we in- 
troduce the Voronoi tessellation as a simple model foam. 
Photon transport in a Voronoi structure using constant 
and thin-film intensity reflectances are discussed in Sees. 
nil Al and IIIIBl respectively. Discussions, conclusions, 
and an outlook are presented in Sec. IIVI 



As a simple model for a three-dimensional disor- 
dered dry foam, we choose the Voronoi tessellation '19']. 
Voronoi foams satisfy the topological requirements on 
edge and face connectivity in Plateau's rules, but not the 
geometric conditions; e.g., the angles between two edges 
are not equal to the tetrahedral angle. To make such a 
Voronoi foam, a distribution of seed points in the sim- 
ulation box is chosen and then Voronoi polyhedrons are 
constructed in complete analogy to the Wigner-Seitz cells 
for periodic lattice sites. We start with a body-centered 
cubic lattice of seed points, which gives the Kelvin foam, 
and then systematically introduce disorder by shifting 
the seed points in random directions. The magnitude 
of the displacement vectors of the seed points is uni- 
formly distributed in the interval [0, h] and disorder in a 
Voronoi foam increases with h. Referring h to the diame- 
ter T^Koivin of a cell in the Kelvin foam [T^ , we study flve 
samples with h/HKc\vin = 0.02,0.05,0.07,0.12, and 0.20. 
Our Voronoi foams are produced by the software Qhull 
|23| . Typically, the simulation box contains around 8200 
cells, 62800 films, and 11200 edges. To simulate the diffu- 
sion of photons in these model foams, periodic boundary 
conditions are implemented. 



In disordered foams several random variables exist; 
e.g., the edge length L, the cell surface S, the cell vol- 
ume V, the number of faces per cell /, and the number of 
edges per face n. A first insight into the cellular structure 
can be gained through the distribution of these variables. 
To characterize foams and their geometrical and topo- 
logical disorder, we have collected data for the following 



quantities: 
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where averages are denoted by ( ) . The first three quan- 
tities give the standard deviations relative to mean val- 
ues of the three respective distributions for L, S, and V. 
The isoperimetric quotient Iq is one for a sphere and, 
therefore, measures how strongly the cells are deformed 
relative to a sphere. Finally, the variances /i2j and fi2,n 
are considered as measures for topological disorder. We 
also note that instead oi AV/V and Iq, alternative mea- 
sures p = (l/)2/3/(y2/3\_ 1 and /3 = (5/(367rl^2)i/3^^ 

are introduced in Ref. (24| . 

It is instructive to compare our samples with slightly 
polydisperse real foams investigated by Monnereau et 
al. [25i] . and random monodisperse foams simulated by 
Kraynik et al. [2g|, see Table HI To produce such foams, 
Kraynik et al. used the Surface Evolver to relax an initial 
Voronoi structure produced from the packing of spheres. 
Another structure that is commonly considered the com- 
plete spatial random pattern, is the Poisson Voronoi tes- 
sellation (PVT) [l^. It allows the analytic evaluation 
of the moments of various distribution functions, and 
consequently has gained much attention. However, PVT 
clearly cannot serve as a model for real foams. As Table 
H] shows, PVT possesses a high degree of topological dis- 
order, ii2j — 11.055 and fi2,n = 2.49, while in real foams 
fi2j < 3.66 and /i2,n < 0.42 
to the simulated foams in [26 1 
ing number of short edges, i.e., the distribution of edge 
lengths reaches its maximum at L = [ig |. 

The distribution of L, 5", V, f and n for our Voronoi 
foams are depicted in Figs. [TJ21 Note that these foams 
are constructed from an initial Kelvin structure, thus 
here Lxoivin, ^Kelvin, and T4:eivin serve as natural units 
of length, surface, and volume, respectively. Characteris- 
tics of the Kelvin foam, our Voronoi foams 1-5, real foams 
[25| , simulated foams [2S,] , and the Poisson Voronoi tessel- 
lation [1^ are summarized in Table HI We shortly discuss 
several points. The standard deviations and variances 
in Table U and the distributions in Figs. [1] and [5] clearly 
demonstrate that the disorder in the Voronoi foams in- 
creases with increasing displacement h of the seed points. 
While Voronoi foam 1 {h / Hy^cWm — 0.02) with its narrow 
distributions is still close to the Kelvin foam, the widths 
of the distributions grow with increasing h and are quite 



25| . Moreover, in contrast 
PVT has an overwhelm- 
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FIG. 1: Distribution of edge length (in units of LkcIv 
the disordered Voronoi foams 1-5. 
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FIG. 2: Distribution of cell surface (in units of S'Koivin), and 
distribution of cell volume (in units of Vkdvin), for the disor- 
dered Voronoi foams 1-5. 



broad in Voronoi foam 5 with h/Hj^civin = 0.20. This 
is especially visible in the distributions for edge length 
L, cell surface S, and cell volume V in Figs. [T] and [2] 
For /i/i?Keivin > 0.20, the Voronoi foam approaches the 
Poisson Voronoi tessellation. The unit cell of a Kelvin 
foam consists of eight hexagonal and six square faces. 
This is especially visible in the bimodal distribution for 
the number of faces with n edges in Fig. [3l which occurs 
in all of our Voronoi foams 1-5. Note that foam 4 ad- 
ditionally contains pentagons and in foam 5 furthermore 
triangles and heptagons appear. In real foams, however, 
the distribution has its single maximum at n = 5. 

Our Voronoi foams constructed from an initial Kelvin 
structure, are indeed distinct from real foams, although 
some of the statistical characteristics of foams 4 and 5 
are close to or agree with the one of real foams. Never- 
theless, as we have just discussed, our Voronoi foams 1-5 
are characterized by increasing topological and geomet- 
rical disorder. In the following, we investigate whether 
this disorder does have any influence on the transport- 
mean-free path I*. 





FIG. 3: Distribution of cells with / faces, and distribution of 
faces with n edges, for the disordered Voronoi foams 1-5. 



III. DIFFUSIVE TRANSPORT OF LIGHT 

We investigate the transport-mean-free path I* for two 
models. 



A. Model of Constant Intensity Reflectance 

To explore the impact of disorder on I*, we assume 
here that liquid films have a constant reflectance r. We 
model single photon paths in a Voronoi foam as a random 
walk with rules motivated by ray optics, i.e., an incom- 
ing light beam is reflected from a face with a probability 
r or it traverses the face with a probability t = 1 — r. 
This naturally leads to a persistent random walk of the 
photons [13i] . where the walker remembers its direction 
from the previous step [23, l2M| . Persistent random walks 
are employed in biological problems [29], turbulent dif- 
fusion [30], polymers 31 [, Landauer diffusion coefficient 
for a one-dimensional solid [32] , and in general transport 
mechanisms [331 ]. More recent applications are reviewed 
in m. 

Our computer program takes 10'* photons at an initial 



position, and launches them in a direction specified by 
polar angles 9 and ip. Then it generates the trajectory of 
each photon following a standard Monte Carlo procedure 
and evaluates the statistics of the photon cloud at times 
T e [10000, 10200,..., 15800] (in units of (L)/c, where c 
denotes the velocity of light). As detailed in Ref. T3|. 
we determine the diffusion constant D from the temporal 
evolution of the average mean-square displacement of the 
photons: (r^) ~ 6Dt. Then the transport-mean-free path 
follows from 



I* 



3D/c. 



(4) 



For angles 6 e [3°, 20°, ..., 71°] and ip G 
[3°, 20°, ..., 156°], the simulation is repeated for each 
intensity reflectance re [0.1, 0.2, ..., 0.9]. As a reasonable 
result, no dependence on the starting point and the 
starting direction is observed. In Fig. |4] we plot the 
transport-mean-free path I* in units of the average cell 
diameter H as a, function of r, where H is defined via 
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FIG. 4: The transport-mean-free path l* (in units of the 
average cell diameter H) as a function of intensity reflectance 
r, for the Kelvin and the disordered Voronoi foams 1-5. The 
Monte Carlo simulation results and the fit r(r)/ H = 0.63(1 — 
r)/r are denoted, respectively, by points and line. 



H = [6(T/)/7r]i/3. The line is a fit to 0.63(1 - r)/r. In 
Fig. m the transport-mean-free paths I* for the regular 
Kelvin foam and the disordered Voronoi foams 1-5 
are not distinguishable from each other. To increase 
the resolution, the rescaled transport-mean-free path 
r(r)/[0.63iJ(l - r)/r\ versus r is plotted in Fig. [3 The 
errorbars reflect the the standard deviations when we 
average over all transport- mean- free paths l*{0,ip) for 
different starting positions and angles. From Fig. [5] 
we find that our numerical results agree well with the 



relation I 



(r) = 0.63iJ-!^(l -bi + h2r) mentioned 



Voronoi 

already in Eq. ([T|) in the introduction. Both constants 
6i and 62 show a slight dependence on disorder in the 
Voronoi foam. The results here already demonstrate 
that disorder in the Voronoi foams does not change I* 
significantly. 



B. Model of Thin-Film Intensity Reflectance 

In the second, more realistic model, we use the inten- 
sity reflectance of thin films. Again, we study the per- 
sistent random walk of photons to probe the influence of 
disorder in the Voronoi foams on the transport-mean-free 
path. 

For a plane wave with wave vector kk incident from 
the air onto a liquid film with normal vector n, the inci- 
dent electric field with a general state of polarization is 
i^incidcnt = ^le^^iei + E2e"f>^e2 + E^e^'t'^e^, where E^ 
and (j)m are respectively the magnitude and phase of the 
field component along the unit vector e^ (m = 1,2,3). 
Taking into account all possible multiple refraction paths 
in the film [35] , the electric field vectors of the transmit- 
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FIG. 5: The transport-mean-free path I* plotted relative to 
0.63H{1 — r)/r as a function of intensity reflectance r. We 
flnd that ivoronoi(»') = 0.6SH{l~r)/r (1 ~ h + b2r) with 
< fei < 0.1 and 62 '^ 0.1. 



ted and reflected waves become 

.^transmitted y^s ''p j (^ .^incident ' OjO -\- ip-tlj\i\c\dcnt 

.C/fcflcctcd — ys r ^p J (^.C/incidcnt ' OjO 7"p£/incidcnt 

+2rp(£;incidont ■ n)n, (5) 

from which one calculates the intensity reflectance r(i) = 
I^EreflcctodP/l^^incidcntP as a function of the incident angle 
i of the plane wave: 



r(i) 



\E, 



incident 



■b\- 



-'-'incident 



iVs 



(6) 



Here | | denotes the magnitude of a complex number, 
and b = h X k. Details of our approach with the explicit 
formulas for the coefficients Tp, tp, r^ and ts are given 
in Ref. [la]. We implemented the reflectance r{i) of Eq. 
^ in our Monte Carlo simulations. Again, we used the 
launching directions for the photons as mentioned in Sec. 
IIII A| but evaluated the statistics of the photon cloud at 
times T £ [30000, 30200,..., 35800] (in units of (L)/c). 
Note that the thin-film reflectance is small except near 
the grazing incidence, thus long simulation times are re- 
quired to achieve the accuracy reported in Sec. IIII Al 

The reflectance r crucially depends on the ratio d/X of 
film thickness d and wavelength A of light, the refractive 
index uq of the film, and the incident angle i. Note that 
even in films as thin as the common black film, r increases 
to 1 close to grazing incidence (i — )■ 90°). This feature 
basically explains why films significantly contribute to 
the diffusion of light in aqueous foams flj, |lfl] . 

First, we assume that all films of the foam have the 
same thickness dav In Fig. [5] we plot I* as a function 
of dav/X for the Kelvin foam and our disordered Voronoi 
foams 1-5. For dav/X < 0.2, the transport-mean-free 
path monotonically decreases as dav increases. There is 
a pronounced minimum around dav/X = 0.2. Between 
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FIG. 6: The transport-mean-free path I* (in units of the av- 
erage cell diameter H) as a function of dav/\ for the Kelvin 
foam and disordered Voronoi foams 1-5. The refractive index 
no of the liquid phase is 1.34. 




FIG. 7: The transport-mean-free path /* of Voronoi foams 1-5 
plotted relative to ^Kcivin ^-s a function of dav/X. For all these 
foams, dw = 0. 



dav/X = 0.2 and 3 oscillations exist that are due to os- 
cillations in the reflectance r of a thin film. The oscilla- 
tions are reduced for more disordered foams but do not 
disappear. The transport-mean-free path I* exhibits two 
strong maxima around dav/X = 0.43 and 0.84. A closer 
inspection of Fig. [5] shows that the heights of the pro- 
nounced maxima at dav/X = 0.43 and 0.84 decrease as 
the disorder increases. To explain this behavior, we note 
that in the Kelvin structure, long straight photon paths 
can cross parallel faces. These long photon paths occur 
for the parallel polarization state since below the Brew- 
ster angle the reflectance is small and especially for the 
most pronounced maximum at dav/X = 0.43 close to zero 
[16j . With increasing disorder in the Voronoi foams these 
parallel faces and, therefore, long straight photon paths. 
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FIG. 8: The transport-mean-free path I* (in units of the aver- 
age ceU diameter H) as a function of dav/X for Voronoi foam 
4 with various thickness distribution. Note, for d^/dav ~ 0.5 
and 0.9, the pronounced maximum disappears. 



disappear. In Fig. [6] the transport-mean-free paths I* 
for the regular Kelvin foam and the disordered Voronoi 
foams 1-5 are not distinguishable from each other in the 
region of dav/^ < 0.3. To increase the resolution, the 
rescaled transport-mean-free path ^*/?Kcivin '^^ Voronoi 
foams 1-5 is plotted for this region in Fig. [71 Clearly, the 
transport-mean-free paths differ from each other by at 
most 10 %. 

Second, in our model we introduce some additional 
randomness in the thickness d of the films assuming that 
it is uniformly distributed in [dav — dw,dav + dw], where 
dav denotes the average thickness and d^ the width of the 
distribution. In Fig. |8]we plot I* as a function of dav/^ 
for Voronoi foam 4 with various thickness distributions. 
Other disordered Voronoi foams show the same behav- 
ior. Obviously, disorder in d decreases the oscillations to 
an approximately constant transport-mean-free path for 
dav/^ > 1- The two strong maxima are reduced notice- 
ably and ultimately disappear for strong disorder in the 
film thickness, as illustrated by the data for d^/dav = 0.5 
and 0.9. However, we still observe that moderate disor- 
der in the thickness {d^/dav ^ 0.2) does not affect I* and 
its monotonic behavior for dav/^ < 0.2. 



IV. DISCUSSION, CONCLUSIONS, AND 
OUTLOOK 



To understand the role of liquid films for light trans- 
port in dry foams, we started with the two-dimensional 
honeycomb [l3| and Voronoi structures [ij]. Then we 
focused on the Kelvin foam [la] to explore the effect of 
space dimension. The Kelvin foam is periodic in space, 
while real foams are disordered. To investigate the in- 



fluence of topological and geometrical disorder on diffu- 
sive photon transport, we have utilized in this paper the 
three-dimensional Voronoi tessellations. 

We have studied the photon's persistent random walk 
in a three-dimensional Voronoi structure based on rules 
motivated by ray optics. In a first model, we used a 
constant intensity reflectance r. The interesting result 
is that the transport-mean-free path for the honeycomb 
structure, the two-dimensional Voronoi foam, the three- 
dimensional Kelvin structure, and the three-dimensional 
Voronoi foam are all determined by the same factor 
(1 — r)/r for constant intensity reflectance r in spite of 
the differences in dimension and structure. For the two- 
dimensional Voronoi foams, we even see the same behav- 
ior as in Eq. ([T]) but with the prefactor 0.63 replaced by 
0.55. Finally, the results here already demonstrate that 
disorder in the Voronoi foams does not change I* signifi- 
cantly [see Figs, m and [5]. This confirms our speculation 
in Refs. [13-.16J that neither the dimension of space nor 
disorder have a strong influence on the magnitude of the 
transport-mean-free path. Note that the factor (1 — r)/r 
expresses the fact that independent of the dimension of 
space and the different shapes of the cells in a foam, pho- 
ton transport is ballistic for r = 0, and that photons stay 
confined to the initial cell for r — 1. 

In a second, more realistic model, we used the intensity 
reflectance of a thin film, with its significant dependence 
on film thickness dav and angle of incidence i. Close 
to grazing incidence {i — >■ 90°), the reflectance always 
sharply increases to one. Thus, a thin film, even as thin 
as the common black film, randomizes the photon direc- 
tion and contributes to I*. Based on extensive Monte- 
Carlo simulations, we paid special attention to the be- 
havior of I* at small dav/X- in real foams dav < 100 nm 
[SS] and, combined with the visible portion of the spec- 
trum (450 nm < A < 750 nm), one arrives at dav/^ < 0.2 
as the relevant region in explaining the measurements of 
the transport- mean- free path in Ref. [9|. Quite inter- 
esting, we found that the monotonic behavior of I* for 
dav/^ < 0.2 does not show a significant dependence on 
the detailed geometrical structure of our Voronoi foams 
and on moderate disorder in the film thickness [see Figs. 
[6118] . 

Real foams obey Plateau's laws. They state that the 
following: i) each film has a constant mean curvature; 
ii) three films meet at angles of 120° at each Pleateau 
border; and iii) these borders meet in fours at the tetra- 
hedral angle [arccos(— 1/3) = 109.47°] to form a vertex. 
A Kelvin foam with fiat faces of its cells approximately 
obeys the second and third law; however, a disordered 
Voronoi foam does not. We found that the transport- 
mean-free path of the Kelvin and the Voronoi structures 
deviate from each other by at most 10% in the experi- 
mentally relevant region for da^/A. This is an indication 
that fulfilling Plateau's second and third laws is of minor 
importance when one tries to understand the role of liq- 
uid films for diffusive transport of light in foams. Even 
the slight curvature of the films might not be very impor- 



tant since in the ray optics approach it just introduces 
some additional randomness in the relevant surface nor- 
mal which varies across the film. Future investigations on 
realistic foam structures, e.g., constructed by the Surface 
Evolver will clarify all these aspects. 

Vera, Saint- Jalmes, and Durian observed the empirical 
law stated in Eq. ^ for 0.008 < e < 0.3. Experimen- 
tal values for I* /H increase from 2 to 20 for decreasing 
e, whereas we determine a range of I* /H between 5 (for 
dav/^ = 0.2) and 25 (for dav/X = 0.06, i.e., for a common 
black film) as illustrated in Fig. [51 That means our the- 
oretical values for I* fall in the same range as the experi- 
mental values for the driest foams in Ref. (9| . This shows 
that liquid films are important for the understanding of 
photon diffusion in dry foams otherwise we should have 
obtained much larger transport-mean-free paths. How- 
ever, to compare our theoretical findings to experiments, 
we should know the contribution of the Plateau borders 
to I*. After all, even in a rather dry foam with a liquid 
fraction of 1 %, most of the liquid resides in these borders 
and their junctions. Furthermore, a relation between the 
film thickness dav and the liquid volume fraction e is 
needed. To the best of our knowledge, systematic mea- 
surements of the film thickness in dry foams have not 
been reported, neither do theoretical models exist. In 
the case of dense oil-in-water emulsions, the dependence 
of the film thickness on the oil- volume fraction has been 



modeled by Buzza and Gates [37|. However their results 
cannot be applied directly to foams. 

In this paper, we focused on the contribution of liquid 



films to the transport-mean-free path of photon diffu- 
sion in dry foams. To achieve a complete understanding 
of the subject, there is still much to do. Of immedi- 
ate interest are modeling of the film thickness variation 
with the liquid volume fraction and following the ran- 
dom walk of photons by taking into account both their 
refraction at thin films and their scattering from Plateau 
borders. To the best of our knowledge, the scattering ma- 
trix for a Plateau border is not available. However, ray 
optics reveals intriguing optical properties of the Plateau 
borders 38, 39]. Using the hybrid lattice gas model for 
two-dimensional foams. Sun and Hutzler [38| emphasize 
that the complex geometry of the Plateau borders affects 
the optical properties of foams. However, Plateau bor- 
ders of a two-dimensional foam do not form a network 
and channeling of the photons, as observed in Ref. |ll| . 
is impossible. Thus an extension of Ref. [38] to three- 
dimensional foams would be interesting. 
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